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[i] Warm pool El Nino (WPEN) events are characterized by positive sea surface 
temperature (SST) anomalies in the central equatorial Pacific. Under present-day climate 
conditions, WPEN events generate poleward propagating wavetrains and enhance 
midlatitude planetary wave activity, weakening the stratospheric polar vortices. The late 
2 1 st century extratropical atmospheric response to WPEN events is investigated using the 
Goddard Earth Observing System Chemistry-Climate Model (GEOSCCM), version 2. 

GEOSCCM simulations are forced by projected late 21st century concentrations of 
greenhouse gases (GHGs) and ozone-depleting substances (ODSs) and by SSTs and sea ice 
concentrations from an existing ocean-atmosphere simulation. Despite known ocean-atmosphere 
model biases, the prescribed SST fields represent a best estimate of the structure of late 21st 
century WPEN events. The future Arctic vortex response is qualitatively similar to that 
observed in recent decades but is weaker in late winter. This response reflects the weaker SST 
forcing in the Nifro 3.4 region and subsequently weaker Northern Hemisphere tropospheric 
teleconnections. The Antarctic stratosphere does not respond to WPEN events in a future 
climate, reflecting a change in tropospheric teleconnections: The meridional wavetrain 
weakens while a more zonal wavetrain originates near Australia. Sensitivity simulations show 
that a strong poleward wavetrain response to WPEN requires a strengthening and 
southeastward extension of the South Pacific Convergence Zone; this feature is not captured 
by the late 21st century modeled SSTs. Expected future increases in GHGs and decreases in 
ODSs do not affect the polar stratospheric responses to WPEN. 

Citation: Hurwitz, M. M., C. I. Garfinkel, P. A. Newman, and L. D. Oman (2013), Sensitivity of the atmospheric response 
to warm pool El Nino events to modeled SSTs and future climate forcings, J. Geophys. Res. Atmos., 118, doi: 10. 1002/ 
2013JD021051. 


1. Introduction 

[ 2 ] Wann pool El Nino (WPEN) events are characterized 
by positive sea surface temperature (SST) anomalies in the 
central equatorial Pacific and have occurred with increasing 
frequency in the past few decades [. Larkin and Harrison, 
2005; Ashok et al., 2007; Kug et al., 2009]. WPEN events 
are distinct from conventional or cold tongue El Nino events, 
typically with stronger, positive SST anomalies in the eastern 
equatorial Pacific [Rasmusson and Carpenter, 1982]. WPEN 
events modulate the extratropical climate [e.g., Ashok et al., 
2007; Garfinkel et al., 2012], 

[ 3 ] In the Southern Hemisphere (SH), WPEN events en- 
hance convective activity in the South Pacific Convergence 
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Zone in austral spring, forcing a tropospheric planetary wave 
that propagates toward SH high latitudes and upward into 
the Antarctic stratosphere [. Hurwitz et al., 2011a, 2011b]. 
This wave enhancement affects Antarctic surface tempera- 
tures [Schneider et al., 2012] and sea ice concentrations [Song 
et al., 2011] and leads to higher polar stratospheric temp- 
eratures and to a weaker polar jet during austral summer, as 
compared with neutral El Nino-Southern Oscillation (ENSO) 
years [Hurwitz et al., 201 la, 201 lb]. Enhanced planetary wave 
driving, higher polar lower stratospheric temperatures, and a 
weaker Antarctic jet have also been identified in response to a 
different definition of WPEN events. Specifically, Zubiaurre 
and Calvo [2012] and Xie et al. [2012] used the El Nino 
Modoki index [Ashok et al., 2007], a different metric by which 
to identify WPEN events, and found SH responses consistent 
with Hurwitz et al. [201 la, 201 lb]. 

[ 4 ] A recent model study by Garfinkel et al. [2012] det- 
ermined that qualitatively, WPEN events have the same 
Northern Hemisphere (NH) impacts as conventional El 
Nino events: A deepened North Pacific tropospheric low 
and enhanced planetary wave driving in boreal winter, lead- 
ing to a weakening of the Arctic vortex. Apparent contradic- 
tions in the findings of Hegyi and Deng [201 1], Xie et al. 
[2012], and Graf and Zanchettin [2012] reflect both limita- 
tions in ENSO analysis when using the short observational 
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Table 1 . Summary of GEOSCCM Simulation Names and Boundary Conditions 3 


Simulation Name 

GHGs 

ODSs 

SSTs and Sea Ice Composites 

Pn 

2005 

2005 

Hadley Centre Global Sea Ice and Sea Surface Temperature version 1 (HadlSSTl) (NTRL) 

P w 

2005 

2005 

HadlSSTl (WPEN) 

IDEALw 

2005 

2005 

HadlSSTl (WPEN in the deep tropical Pacific; NTRL elsewhere) 

Fn 

2100 

2100 

Community Climate System Model version 3.0 (CCSM3) 2070-2100 (NTRL) 

F w 

2100 

2100 

CCSM3 2070-2100 (WPEN) 

FSST n 

2005 

2005 

CCSM3 2070-2100 (NTRL) 

FSST W 

2005 

2005 

CCSM3 2070-2100 (WPEN) 

fcli n 

2100 

2100 

HadlSSTl (NTRL) 

FCLIw 

2100 

2100 

HadlSSTl (WPEN) 

A n 

2005 

2100 

HadlSSTl (NTRL) 

A w 

2005 

2100 

HadlSSTl (WPEN) 

Bn 

2100 

2100 

Hadley Centre Global Environmental Model, version 1 (HadGEMl) 2070-2100 (NTRL) 

B w 

2100 

2100 

HadGEMl 2070-2100 (WPEN) 

C N 

2005 

2005 

CCSM3 2000-2030 (NTRL) 

Cw 

2005 

2005 

CCSM3 2000-2030 (WPEN) 


“All are 50 year time slice simulations, as described in the text. 


record and the sensitivity of the NH response to the precise 
definition of central Pacific El Nino events [ Garfinkel 
et al., 2012]. 

[ 5 ] Dramatic polar stratospheric climate change is exp- 
ected by the end of the 21st century. Following moderate 
emissions scenarios, the atmospheric C0 2 concentration will 
approximately double while concentrations of ozone-depleting 
substances (ODSs) will be greatly reduced. While greenhouse 
gas (GHG)-related climate change will moderately cool 
the polar lower stratosphere in autumn, decreasing levels 
of ODSs will lead to ozone recovery and thus to warming 
trends in spring [Hitchcock et al., 2009; Hurwitz and 
Newman, 2010], Late 21st century changes in El Nino 
amplitude and frequency are unclear: The Coupled Model 
Intercomparison Project (CMIP5) models suggest no future 
trend in the overall ENSO amplitude [Stevenson, 2012], 
while other model studies predict that the pattern of 2 1 st 
century SST trends will favor central Pacific wanning 
[Yeh et al., 2009; Xie et al., 2010]. If so, in the polar strato- 
sphere, the dynamical wanning associated with frequent 
WPEN events would offset some of the predicted radiative 
cooling by GHGs. 

[6] Not only may the frequency of WPEN events 
change in a future climate but also may the tropospheric 
planetary wave response. Several model studies predict that 
the NH circulation response (i.e., teleconnections) to strong 
ENSO events will change under enhanced CO 2 conditions 
[e.g., Collins, 2000; Muller and Roeckner, 2006; Kug 
et al., 2010]. 

[ 7 ] Since WPEN events have contributed to polar strato- 
spheric variability in recent decades, a thorough evaluation 
of the extratropical impacts of these events in a late 21st 
century climate will improve overall estimates of future 
variability. This paper describes the future polar stratospheric 
response to WPEN events, as simulated by the Goddard Earth 
Observing System Chemistry-Climate Model (GEOSCCM), 
and uses idealized and sensitivity simulations to pinpoint 
the causes of potential changes in the WPEN response 
due to modeled SST fields and future climate conditions. 
The model and simulations are described in section 2. 
Section 3 presents the results, comparing WPEN and 
ENSO neutral simulations, under present, future, and ideal- 
ized climate conditions. Section 4 presents a brief summary 
and discussion. 


2. Model Description and GEOSCCM 
Simulations 

2.1. Model Description 

[8] The Goddard Earth Observing System Chemistry-Climate 
Model, (GEOSCCM) version 2 couples the GEOS-5 atmo- 
spheric general circulation model (GCM) [Rienecker et al., 
2008; Molod et al., 2012] and a comprehensive strato- 
spheric chemistry scheme [Kawa et al., 2003; Pawson 
et al., 2008], The model is run at 2° latitude x 2.5° longitude 
horizontal resolution and 72 vertical layers, with a model 
top at 0.01 hPa. Predicted distributions of water vapor, ozone, 
the primary greenhouse gases (C0 2 , CH 4 , and N 2 0), CFC- 
11, CFC-12, and HCFC-22 feed back to the radiative calcula- 
tions. The GEOSCCM performed well in the Stratospheric 
Processes and their Role in Climate (SPARC) Chemistry- 
Climate Model Validation (CCMVal) [2010] detailed evalua- 
tion of stratospheric processes. 

[ 9 ] The present study uses the fonnulation of the GEOSCCM 
as documented by Hurwitz et al. [2011b]. Specifically, the 
stratospheric chemistry module is the same as that documented 
by SPARC CCMVal [2010], while the GEOS-5 GCM has been 
updated. The updated GCM generates a spontaneous quasi- 
biennial oscillation, with a realistic period (27 months at 
30hPa and 50hPa [Hurwitz et al., 2013]) and zonal wind 
amplitude, and improved tropospheric stationary wave pat- 
terns in the Southern Hemisphere. This formulation of the 
GEOSCCM has been successfully used to evaluate the re- 
sponse of the Antarctic stratosphere to warm pool El Nino 
events in a present-day climate [Hurwitz et al., 2011b], 
compare the NH responses to central Pacific and eastern 
Pacific El Nino events [Garfinkel et al., 2012], and isolate 
the impact of North Pacific SSTs on the Arctic winter cli- 
mate [Hurwitz et al., 2012], 

2.2. GEOSCCM Simulations 

[ 10 ] Nine GEOSCCM simulations, plus six supplementary 
simulations, are used to isolate the impact of WPEN events 
on the upper troposphere and polar stratosphere, in present- 
day and likely future climate conditions (Table 1). All simu- 
lations are 50-year time slices, with fixed concentrations of 
the primary greenhouse gases (GHGs) and ozone-depleting 
substances (ODSs), and sea surface temperature (SST) and 
sea ice fields with repeating annual cycles. This experimental 
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Figure 1 . F N -P N global and annual mean temperature 
differences (K). 

design provides large samples of the atmospheric response to 
each set of climate conditions. Pairs of simulations each con- 
tain a WPEN simulation and an ENSO neutral (NTRL) sim- 
ulation, with SST and sea ice fields generated by compositing 
WPEN or NTRL events, respectively. Each ENSO event 
included in the climatologies spans from the July preceding 
the boreal winter peak in tropical SST anomalies through 
June of the following year. 

[ 1 1 ] A first pair of simulations represents the present-day 
climate (“P”). As reported by Hurwitz et al. [2011b], 2005 
GHG and ODS concentrations are prescribed in both simu- 
lations, but the prescribed SST and sea ice conditions are 
distinct. The 1991-1992 and 1994-1995 WPEN events are 
composited to create the SST and sea ice climatologies used 
as boundary conditions for the P w simulation; SST and sea 
ice concentrations are derived from HadlSSTl [Rayner et al., 
2003]. Similarly, 10 NTRL events spanning the 1979-2011 
period are composited to create SST and sea ice climatologies 
used as boundary conditions in the P N simulation. 

[ 12 ] An additional, idealized WPEN simulation (“IDEALw”; 
“CPWideal” in Garfinkel et al. [2012]) is compared with P N . 
Differences between 1DEAL W and P N isolate the contribution 
of deep tropical central Pacific SST anomalies to the atmo- 
spheric WPEN teleconnections. In IDEAL W , SSTs between 
10°N, 10°S, 140°E, and 240°E are identical to those in the 
P w simulation. SSTs poleward of 20° latitude, west of 115°E 
and east of 265°E are identical to those in the P N simulation. 
Between these two regions, SSTs in 1DEAL W are linearly 
interpolated between the P w and P N climatologies. 

[ 13 ] A second pair of simulations represents the expected 
late 21st century or “future” climate (“F”). Concentrations 
of GHGs [A1B scenario; Solomon et al., 2007] and ODSs 
[Al scenario; World Meteorological Organization, 2011] 


are both fixed at 2100 values. SSTs and sea ice are taken 
from the Community Climate System Model, version 3.0 
(CCSM3) (http://www.ccsm.ucar.edu), of the Coupled 
Model Intercomparison Project, phase 3 (CM1P3) (http:// 
www-pcmdi.llnl.gov/ipcc/about_ipcc.php) A1B scenario 
simulations of the 21st century, for the 2070-2100 period. 
WPEN and NTRL events are identified using similar criteria 
as for observed events: WPEN events are identified when the 
September-October-November (SON) seasonal mean Nino 4 
index (monthly SST anomalies in the 5°S-5°N, 160°-210°E 
region; http://www.cpc.ncep.noaa.gov/data/indices/) ex- 
ceeds one standard deviation from the 2070-2100 mean 
and is larger than the equivalent Nino 3 index (monthly 
SST anomalies in the 5°S-5°N, 210°-270°E region). 
NTRL events are identified when the SON seasonal mean 
Nino 3 and Nino 4 indices are both within 0.7 standard 
deviations from the respective 2070-2100 means. The 6 
WPEN and 12 NTRL events are composited to create SST 
and sea ice boundary conditions, respectively, for the F w 
and F n simulations. 

[ 14 ] The future simulations represent large changes in the 
global and annual mean temperature from the present day. 
Specifically, Figure 1 shows differences between the NTRL fu- 
ture simulation and the NTRL present simulation (i.e., F N -P N ). 
The troposphere warms by approximately 2 K, while strong 
stratospheric cooling increases with height, exceeding 7 K in 
the upper stratosphere. 

[ 15 ] Differences in tropical Pacific SSTs, as prescribed in 
the WPEN and NTRL simulations, are shown in Figure 2. 
By definition, WPEN-NTRL SST differences are positive 
in the central equatorial Pacific (i.e., Nino 4 region), in the 
present-day (Figures 2a and 2d), idealized present-day 
(Figures 2b and 2e), and future (Figures 2c and 2f) cases, 
with SST differences exceeding 1 K in some regions in 
boreal winter. SON (December-January-February (DJF)) 
seasonal mean Nino 4 anomalies are comparable in the three 
sets of SST anomalies: 0.76 K (0.89 K) in the P simulations, 
0.73 K (0.83 K) in IDEAL W , and 0.76 K (0.83 K) in the F 
simulations (Table 4). Outside of the deep tropical Pacific, 
SST anomalies differ between the various simulations. 

[ 16 ] Two additional pairs of simulations separate the im- 
pacts of changing radiative forcing and modeled late 21st 
century SSTs on the polar stratospheric response to WPEN 
(Table 1). In the “FCLI” simulations, present-day SSTs are 
prescribed (Figures 2a and 2d), with 2100 values for the 
GHGs and ODSs. In the “FSST” simulations, modeled late 
21st century SSTs are prescribed (Figures 2c and 2f), with 
2005 values for the GHGs and ODSs. 

[ 17 ] Six supplementary simulations test the sensitivity of 
the WPEN teleconnections to ozone recovery (section 3.3) 
and to the choice of prescribed SST fields (section 3.4). 
The “A” pair of simulations (Table 1) tests the sensitivity to 
ozone recovery. These simulations are forced with observed 
SSTs and 2005 values of GHGs, but with projected 2100 
values of ODSs. The “B” pair of simulations (Table 1) is 
identical to the F simulations, except that the late 2 1 st century 
SSTs are taken from a Hadley Centre Global Environmental 
Model, version 1 (HadGEMl) simulation, as is further 
discussed in section 3.4. The “C” pair of simulations is 
identical to the P simulations, except that early 21st century 
modeled SSTs from CCSM3 are used in place of 
observed SSTs. 
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Figure 2. WPEN-NTRL SST differences in the (a-c) SON season and (d-f) DJF season. SST differences 
prescribed in the (a and d) P and FCLI simulations, (b and e) IDEALw - Pw, and (c and f) F and FSST 
simulations. 


3. Results 

3.1. Polar Stratospheric Response to WPEN 

[is] In the present-day simulations, the polar vortex 
weakens in response to WPEN events, i.e., positive WPEN- 
NTRL geopotential height differences poleward of 60° latitude 
in both the Arctic (Figure 3a) and Antarctic (Figure 3b). The 
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Arctic response to WPEN events is strongest in late boreal 
winter [Garfinkel et al., 2012], while in the Antarctic the 
response to WPEN events is strongest in austral spring 
[Hurwitz et al., 201 lb]. In the upper stratosphere, the magni- 
tude of the Arctic response is greater than the Antarctic 
response. In the Arctic, these geopotential height differences 
correspond with lower stratospheric polar cap temperature 


Antarctic 
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Month 



Figure 3. WPEN-NTRL polar cap geopotential height differences, in October through March, as a func- 
tion of altitude (m). (a and b) P simulations, (c and d) F simulations, (e and f) F-P differences. Note the 
different color scales for the Arctic in Figures 3a, 3c, and 3e and Antarctic in Figures 3b, 3d, and 3f. 
White contours indicate zero difference. Black crosses indicate differences significant at the 95% level, 
in a two-tailed t test. 
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Figure 4. WPEN-NTRL polar cap temperature differences, in October through March, as a function of 
altitude (K). (a and b) P simulations, (c and d) F simulations. White contours indicate zero difference. 
Black crosses indicate differences significant at the 95% level, in a two-tailed t test. 


increases of approximately 1 K in the early and midwinter, 
and a larger 3-4 K wanning in February and March 
(Figure 4a). The Antarctic lower stratosphere warms by 
1-2 K, with the largest temperature differences in October 
through December (Figure 4b). 

[ 1 9 ] The polar stratospheric response to WPEN events changes 
in a future climate. Figures 3c, 3d, 4c, and 4d show the equiv- 
alent WPEN responses in the future simulations, while 
Figures 3e and 3f show differences between the future and 
present-day geopotential height responses. In the Arctic, 
future WPEN events continue to weaken and warm the polar 
vortex (Figures 3c and 4c), though the response is weaker than 
that simulated in a present-day climate in October, February, 
and March (Figure 3e). In the lower stratosphere, the polar 
cap temperature response to WPEN is approximately 1-2 K. 
In the F simulations, warming of the Antarctic lower strato- 
spheric response is replaced by weak cooling in October and 
a null response in later months (Figures 3d and 4d). The 
remainder of section 3 will make use of the idealized and sen- 
sitivity simulations to explain the simulated future changes in 
the polar stratospheric responses to WPEN. 

3.2. Tropospheric WPEN Teleconnections 

[ 20 ] Future changes in the polar stratospheric response to 
WPEN events reflect changes in tropospheric planetary wave 
driving, which in turn reflect changes in the tropospheric 
circulation response pattern (or “teleconnections”) to WPEN 
events. Because the Arctic stratospheric response to WPEN 
is strongest in late boreal winter, the Northern Hemisphere 
(NH) tropospheric responses are examined for the preceding 
months (i.e., the DJF season). Because the Antarctic strato- 
spheric response is strongest in austral early summer, the 
Southern Hemisphere (SH) tropospheric responses are exam- 
ined for the SON season. 

[ 21 ] Under present-day conditions, WPEN events generate 
a strong poleward, upper tropospheric wavetrain response 


in the NH in boreal winter, i.e., increased geopotential 
heights in the tropical central Pacific, a low centered at 40°N, 
220°E, and increased heights centered at 60°N, 250°E 
(Figure 5a; see also Garfinkel et al. [2012]). The first tempo- 
ral empirical orthogonal function (EOF) of WPEN-NTRL 
differences in the P simulations highlights this poleward 
wavetrain pattern (Figure 5e). Planetary wave activity in the 
North Pacific is enhanced during observed WPEN events. 
Figure 6 shows the results of Rossby wave source calcula- 
tions [based on Jin and Hoskins , 1995] derived from the 
simulated upper tropospheric wind fields. The average of 
the 250 hPa and 300 hPa pressure levels is shown [as in 
Hurwitz et al., 2011b], as the 250 hPa and 300 hPa levels 
emphasize differences at mid-latitudes. Rossby wave source 
is a concise way of determining the forcing by Rossby waves, 
due to the divergent component of the flow. The deepened 
North Pacific low is located downstream of a Rossby wave 
source (Figure 6a). Another measure of tropospheric plane- 
tary wave driving is the December/January mean eddy heat 
flux at 40-80°N, lOOhPa [Newman et al., 2001]: Eddy heat 
flux is enhanced in response to WPEN events (i.e., P w - Pn), 
in a contemporary climate (Table 2). 

[ 22 ] Under future climate conditions, the NH tropospheric 
response retains the same pattern but weakens, consistent 
with the polar stratospheric response. Specifically, in the 
F simulations, the upper tropospheric geopotential height 
response to WPEN strongly resembles the response in the P 
simulations (compare Figures 5a and 5d). The first EOF of 
the F w — F n response (Figure 5f) is indistinguishable from 
the first EOF of P w — P n (Figure 5e); the two patterns are 
correlated with r = 0.92. In the future simulations, WPEN 
enhances tropospheric wave driving (Table 2); however, the 
North Pacific Rossby wave source response to WPEN is 
weaker in the future simulations (Figure 6d), as are the 
geopotential height differences at 60°N, 250°E, than in the 
present-day simulations. 
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Figure 5. (a-d) WPEN-NTRL 250 hPa geopotential height differences (m) in the DJF season in the (a) P 
simulations, (b) FSST simulations, (c) FCLI simulations, and (d) F simulations. White contours indicate zero dif- 
ference. Black crosses indicate differences significant at the 95% level, in a two-tailed t test, (e-f) Temporal first 
EOFs of WPEN-NTRL geopotential height differences in the North Pacific, at 250 hPa: P and F in Figures 5e 
and 5f, respectively; thick, black arrows indicate the approximate direction of the poleward wavetrains. 
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Figure 6. WPEN-NTRL differences in DJF upper tropospheric Rossby wave source [10 10 s -2 ], as described in 
the text, in (a) P simulations, (b) FSST-P differences, (c) FCL1-P differences, and (d) F-P differences. Black crosses 
indicate differences significant at the 95% level, in a two-tailed t test. White contours indicate zero difference. 
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Table 2. Magnitude of the December/January Eddy Heat Flux at 
40-80°N, 100 hPa a 


Pair of Simulations 

WPEN 

NTRL 

p 

11.01 ±0.22 

10.28±0.25 

F 

14.01 ±0.38 

13.21 ±0.35 

FCLI 

13.50±0.40 

1 1 .88 ± 0.34 

FSST 

14.43±0.37 

13.42±0.45 


a WPEN-NTRL differences are statistically significant at the 95% level in 
the P, FCLI, and FSST simulations, and at the 85% level in the F simulations, 
based on two-tailed t tests. 


[ 23 ] Sensitivity simulations suggest that it is the change in 
the prescribed SST field, rather than the change in ODS or 
GHG concentrations that modifies the stratospheric response 
to WPEN events. In the FCLI simulations, where present-day 
SST climatologies are prescribed under future GFIG and 
ODS concentrations, Rossby wave source differences are in- 
distinguishable from the present-day differences (Figure 6c). 
In the FSST simulations, where future SST climatologies are 
prescribed under present-day GFIG and ODS concentrations, 
the North Pacific Rossby wave source weakens (Figure 6b 
and Table 2), as in the F simulations. 

[ 24 ] Weaker upper tropospheric and Arctic stratospheric 
responses to WPEN, in the simulations where late 2 1 st cen- 
tury modeled SSTs are prescribed, are caused by weaker 
tropical SST forcing and not by future SST changes and/or 
model biases outside of the tropical Pacific. Table 4 shows 
that while the prescribed SST differences in the Nino 4 
region in various simulations are comparable for the D JF sea- 
son (i.e., 0.89 K and 0.83 K), as chosen by the experimental 
design, SST differences directly to the east, i.e., in the Nino 
3.4 region (5°S-5°N, 190°-240°E; http://www.cpc.ncep. 
noaa.gov/data/indices/) are nearly twice as large in the simu- 
lations forced by observed SSTs (i.e., 1.60 K) than in the 
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simulations forced by modeled late 21st century SSTs 
(0.90 K). That is, the magnitude of the extratropical atmo- 
spheric response to WPEN depends on the SST forcing in the 
Nino 3.4 region. Conversely, the WPEN response is not sensi- 
tive to extratropical SST differences: The present-day structure 
and magnitude of the upper tropospheric response and subse- 
quent warming of the Arctic stratosphere (Figures 4a and 5a) 
are seen in the idealized WPEN simulation (Figure 7). 

[ 25 ] As for the NH response, WPEN events generate a 
strong poleward, upper tropospheric wavetrain response in 
the South Pacific in austral spring. Geopotential height 
increases in the tropical, central South Pacific, decreases 
around 40°S, 210°E, and increases around 60°S, 240°E 
(Figure 8a; see also Hurwitz et al. [201 la, 201 lb]). The first 
temporal empirical orthogonal function (EOF)) of WPEN- 
NTRL differences in the P simulations highlights this pole- 
ward wavetrain pattern (Figure 8e). This pattern resembles 
the Pacific-South America- 1 (PSA-1) pattern, as identified 
by Mo [2000] and Vera et al. [2004], and is the second lead- 
ing mode of observed interannual variability in the South 
Pacific region. As for the North Pacific low, the deepened 
South Pacific low at 40°S, 210°E is fed by an anomalous 
Rossby wave source (Figure 9a). The Rossby wave source 
results from the strengthening of convective activity in the 
southeastern South Pacific Convergence Zone (SPCZ), in re- 
sponse to WPEN events [Hurwitz et al., 2011a], Figure 10a 
shows the WPEN-related decrease in outgoing longwave 
radiation (OLR) (a proxy for increased convection) in the trop- 
ical central Pacific and in the diagonal region between 15°S, 
220°E, and 40°S, 270°E, in the P simulations. Furthermore, 
at high latitudes, the structure of the geopotential height anom- 
alies is in phase with the climatological field [Zubiaurre and 
Calvo, 2012] enhancing wave driving to the stratosphere: 
October/November eddy heat flux at 40-80°S, lOOhPa 
increases in response to WPEN (Table 3). 


Southern Hemisphere 
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0 60 120 180 240 300 360 
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Figure 7. (a and b) IDEAL W -P N differences in polar cap temperature, in October through March, as a 
function of altitude (K). (c) DJF geopotential height anomalies at 250 hPa and (d) SON geopotential height 
anomalies at 250 hPa. White contours indicate zero difference. Black crosses indicate differences signifi- 
cant at the 95% level, in a two-tailed t test. 
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Figure 8. (a-d) WPEN-NTRL 250 hPa geopotential height differences (m) in the SON season in the (a) P 
simulations, (b) FSST simulations, (c) FCLI simulations, and (d) F simulations. White contours indicate zero 
difference. Black crosses indicate differences significant at the 95% level, in a two-tailed t test, (e-f) Temporal 
first EOFs of WPEN-NTRL geopotential height differences in the South Pacific, at 250 hPa: P, with dominant 
poleward wavetrain (PSA-1) in Figure 8e and F with dominant zonal wavetrain (PSA-2) in Figure 8f. Thick, 
black arrows indicate the approximate direction of the PSA- 1 -like and PSA-2 -like wavetrains, respectively. 


[ 26 ] When late 21st century SSTs and 2 100 climate condi- 
tions are prescribed (Figure 8d), the structure of the pole- 
ward wavetrain response remains but with significantly 
smaller geopotential height differences at 40°S and 60°S. 
In addition, a region of positive height differences is seen 
around Australia (30°S, 120°E), part of a more zonal 
wavetrain with negative height differences around 60°S, 
150°E, positive height differences around 60°S, 240°E, 
and negative differences around 70°S, 300°E. The dominant 
mode of variability of the Fw— F n differences i.e., the 1st 
temporal empirical orthogonal function (EOF)) highlights 
this zonal wavetrain pattern (Figure 8f). The 1st EOF pat- 
tern resembles the Pacific-South America-2 (PSA-2) mode 
of observed variability [Mo, 2000; Vera et al., 2004], 
Corresponding with the changes in the upper tropospheric 
teleconnection pattern, the Rossby wave source at 30°S, 
240°E weakens in the F simulations relative to the P simula- 
tions, and a relative Rossby wave source appears around 30° 
S, 11 0°E (Figure 9d). 

[ 27 ] As for the Northern Hemisphere teleconnections, 
sensitivity simulations suggest that it is the change in South 
Pacific SSTs rather than the change in ODS or GHG 
concentrations that modifies the tropospheric response and 


ultimately the stratospheric response to WPEN events in the 
SH. In the FCLI simulations, the WPEN teleconnection pattern 
resembles the present-day response (compare Figures 8a and 
8c), with a strong poleward (i.e., PSA-1 -like) geopotential height 
pattern. The pattern of the temporal first EOF of WPEN-NTRL 
differences in the FCLI simulations is strongly correlated with 
that in the P simulations (r = 0.94). As in the present-day 
simulations [Hurwitz et al., 201 lb], strong poleward planetary 
wave activity leads to a significant increase in eddy heat flux 
at 40-80°S, lOOhPa in October/November during WPEN 
events (Table 3). The South Pacific Rossby wave source at 
30°S, 240°E in fact strengthens slightly (Figure 9). 

[ 2 s] The tropospheric response to WPEN changes when 
late 21st century modeled SSTs are prescribed. In the FSST 
simulations, the WPEN teleconnection pattern most resembles 
the future climate response: a weaker PSA-l-like response 
as well as a strong, zonal PSA-2-like response (compare 
Figure 8). The pattern of the temporal first EOF in the FSST 
simulations is strongly correlated with that in the F simulations 
(r=0.91). Weakening ofthe PSA-l-like response weakens the 
tropospheric planetary wave driving of the stratosphere, lead- 
ing to insignificant eddy heat flux responses to WPEN events 
in the FSST and F simulations (Table 2). 
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Figure 9. WPEN-NTRL differences in SON upper tropospheric Rossby wave source [10 10 s -2 ], as described in 
the text, in (a) P simulations, (b) FSST-P differences, (c) FCL1-P differences, and (d) F-P differences. Black crosses 
indicate differences significant at the 95% level, in a two-tailed t test. White contours indicate zero difference. 


[ 29 ] A strong meridional wavetrain response requires posi- 
tive SST forcing from the south-central Pacific, in the region 
200-240°E, 10-20°S. The present-day and idealized simu- 
lations are forced by positive SST anomalies in the south- 
central Pacific (180-220°E, 5-20°S; see Figures 2a and 2b), 
while SST anomalies in the simulations forced by modeled late 
21st century SSTs lack sufficiently positive SST anomalies in 
this region (Figure 2c). This difference in SST forcing leads to 
a weaker convective response in the southeastern South 
Pacific Convergence Zone (i.e., the region of 240°E, 30°S) than 
in the present-day simulations (compare Figures 10c and lOd 
with Figures 10a and 10b), and thus to a much weaker Rossby 
wave source (e.g., Figures 9b and 9d), weaker meridional 


wavetrain response (Figures 8b and 8d) and ultimately to the 
failure to wann the Antarctic lower stratosphere (compare 
Figure 4d to Figures 4b and 7). 

[ 30 ] The appearance of the zonal PSA-2-like wavetrain re- 
sponse in the simulations forced by late 21st century SSTs 
(see Figure 8), with positive geopotential height differences 
at 30°S, 120°E, suggests a heating anomaly near Australia. 
Figure 1 1 shows the 2 m (i.e., surface) temperature responses 
to WPEN events in the present-day, idealized and future sim- 
ulations in austral spring. Note the striking Australian tem- 
perature anomaly exceeding 2K in the F simulations, as 
compared with the present-day and idealized simulations. 
These Australian heating anomalies are seen throughout the 




150 180 210 240 270 300 

Longitude 



150 180 210 240 270 300 

Longitude 


Figure 10. WPEN-NTRL differences in SON seasonal mean OLR (W m“ 2 ) for (a) P W -P N , (b) IDEAL W - 
P N , (c) FSST w -FSST n , and (d) F W -F N . Black crosses indicate differences significant at the 95% level, in a 
two-tailed t test. White contours denote zero difference. 
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Table 3. Magnitude of the October/November Eddy Heat Flux at 
40-80°S, 100 hPa a 


Pair of Simulations 

WPEN 

NTRL 

p 

7.32 ±0.22 

6.80±0.26 

F 

7.08±0.19 

6.89 ±0.24 

FCLI 

7.26 ±0.22 

6.39 ±0.22 

FSST 

6.99 ±0.20 

6.86 ±0.20 


a WPEN-NTRL differences are statistically significant at the 85% (99%) 
level in the P (FCLI) simulations, based on two-tailed 1 tests. 


tropospheric column, i.e., consistent with the appearance of a 
Rossby wave source at 30°S, 1 10°E (Figures 9b and 9d). 

3.3. Sensitivity to Ozone Recovery 

[ 31 ] The previous section showed that in the GEOSCCM, 
the extratropical atmospheric response to WPEN is insensi- 
tive to predicted future changes in GHG and ODS concentra- 
tions. The WPEN response might lack sensitivity to climate 
change because the impacts of simultaneously decreasing 
ODSs (ozone recovery, weakening of the polar vortices and 
equatorward shifting of the tropospheric midlatitude jet) 
and increasing GFIGs (strengthening of the polar vortices 
and poleward shifting of the tropospheric jet) may cancel 
out each other [Polvani et al., 2011; Arblaster et al., 2011], 
To test this hypothesis, an additional pair of GEOSCCM 
simulations (“A”; Table 1) was performed. This pair of sim- 
ulations isolates the impact of ozone recovery on the polar 
stratospheric response to WPEN events, by prescribing 
present-day SSTs and GFIGs but projected 2100 values for 
the ODSs. 

[ 32 ] In isolation, ozone recovery has a negligible impact on 
the tropospheric and stratospheric responses to WPEN (not 
shown). While Antarctic ozone recovery (i.e., additional ra- 
diative heating by stratospheric ozone) warms the polar 
stratosphere, it has no impact on the Antarctic stratospheric 
response to WPEN : The response under ozone recovery is in- 
distinguishable from the present-day simulations. 

3.4. Sensitivity to the Prescribed SST Fields 

[ 33 ] The results presented in sections 3.1 and 3.2 suggest 
that the atmospheric response to WPEN is sensitive to the 
prescribed SST field. ITowever, only one pair of modeled 
SST fields was tested in the set of simulations as described 
above. Two additional pairs of simulations were performed 
for comparison with the F simulations. Together, these three 
pairs of simulations measure, to some degree, the sensitivity 
of the results to the prescribed set of modeled SSTs. 

[ 34 ] The first pair of simulations (“B”; Table 1) is forced 
with 2100 values of GFIGs and ODSs, and late 21st century 


SSTs. SSTs and sea ice are taken from a Hadley Centre 
Global Environmental Model, version 1 (HadGEMl) simula- 
tion [Johns, 2006] using the 2070-2100 period of a A1B 
scenario simulation. WPEN and NTRL events are identified 
as for the F simulations; 6 WPEN and 10 NTRL events are 
composited to create the SST and sea ice boundary con- 
ditions. WPEN— NTRL SST anomalies in the Nino 4 and 
Nino 3.4 regions are comparable to those in the F simulations 
(Table 4). Note that future simulations forced with HadGEMl 
SSTs [Johns, 2006] yield nearly identical changes in global 
mean temperature, despite the higher climate sensitivity to 
doubled C0 2 in HadGEMl than in CCSM3 [Kiehl et al., 
2006], due to cold biases in the present-day climate. 

[ 35 ] The second pair of simulations (“C”; Table 1) is 
forced with 2005 values of GHGs and ODSs, and early 
21st century SSTs. SSTs and sea ice values are taken from 
the CCSM3 A1B scenario simulation, as in the F simulations, 
but for the 2000-2030 period. The 4 WPEN and 10 NTRL 
events are composited to create the SST and sea ice 
boundary conditions. 

[36] The B and C pairs of simulations, forced with 
HadGEMl and CCSM3 SSTs, yield qualitatively similar 
atmospheric responses to WPEN as in the F simulations. 
The poleward, upper tropospheric wavetrains are weaker 
than in the P simulations, weakening planetary wave driving 
of the stratosphere. Compared with the P simulations, the 
Arctic vortex warms less in late boreal winter. This result is 
consistent with the weaker SST anomalies in the Nino 3.4 
region during the DJF season (0.57 K and 0.61 K, in the two 
additional pairs of simulations, respectively) as compared 
with the observed anomalies (1.60 K) (Table 4). Figure 12a 
shows the positive relationship between SST anomalies in 
the Nino 3.4 region and the late winter Arctic stratospheric 
temperature response. 

[ 37 ] The Antarctic stratosphere does not warm in the B and 
C simulations with modeled SSTs. In the C simulations, the 
magnitude of the SST anomalies in the south-central Pacific 
is weaker than observed, leading to a weak convective 
response in the SPCZ region and consequently to a weak plan- 
etary wave response at high latitudes (Figure 12b). In the B 
simulations, the magnitude of the south-central Pacific SST 
anomaly resembles that in the observations. However, in B N , 
the high-latitude stationary wave patterns are not well simu- 
lated. Destructive interference between the patterns in B N 
and B w inhibits the propagation of the wave response to high 
latitudes and ultimately to the Antarctic stratosphere (i.e., the 
weak polar stratospheric temperature response indicated by 
the white square in Figure 12b). A zonal PSA-2-like wavetrain 
does appear in the B simulations, similarly to that seen in the F 
simulations, but with a smaller amplitude. 



Figure 11. WPEN-NTRL 2 m temperature differences (K) in the equatorial and South Pacific in austral 
spring, for (a) P W -P N , (b) IDEAL W -P N , and (c) F W -F N . Black crosses indicate differences significant at 
the 95% level, in a two-tailed t test. White contours denote zero difference. 
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Table 4. SST Anomalies in the Nino 4 and Nino 3.4 Region 


Pairs of Simulations 

SON 

DJF 

Niiio^l 

Nino-3.4 

Nino— 4 

Nino-3 .4 

P, FCLI, A 

0.76 

0.97 

0.89 

1.60 

IDEALyy-PiM 

0.73 

0.96 

0.83 

1.54 

F, FSST 

0.76 

0.90 

0.83 

0.90 

B 

0.96 

0.89 

0.73 

0.57 

C 

0.73 

0.92 

0.69 

0.61 


4. Summary and Discussion 

[38] The late 21st century atmospheric response to warm 
pool El Nino events was investigated with a set of 
GEOSCCM simulations. Pairs of 50 year time slice provided 
large samples of the response to ENSO neutral and WPEN 
conditions, in both present-day and in projected late 21st 
century climate conditions. The GEOSCCM can simulate 
the observed, present-day response to WPEN events when 
forced by observed SST fields [Hurwitz et al., 2011b; 
Garfinkel et al., 2012]. For the GEOSCCM simulations of 
the future, late 21st century sea surface temperatures and 
sea ice concentrations, an existing ocean-atmosphere simula- 
tion were prescribed as boundary conditions (Figure 2), in 
addition to projected 2100 concentrations of the primary 
greenhouse gases and ozone-depleting substances. Despite 
known ocean-atmosphere model biases, the prescribed SST 
fields represent a best estimate of the structure of late 2 1 st 
century WPEN events. In the Arctic, the polar vortex weak- 
ened in response to WPEN events, as for conventional El 
Nino events [Herceg Bulic et al., 2012], though the late win- 
ter response to WPEN events was not as strong as that in a 
contemporary climate. The Antarctic vortex weakened in 
response to WPEN events in the present-day climate but 
did not weaken or warm in response to WPEN events in a 
future climate (Figures 3 and 4). 

[ 39 ] The combined impact of increasing GFIGs and de- 
creasing ODSs on the polar stratosphere on the response to 
WPEN was negligible, despite large differences between 
the present-day and projected future climates (Figure 1). 
That is, when present-day SSTs were prescribed, the polar 
response to WPEN (both in the Arctic and Antarctic) was 
the same in both present-day and future climate conditions: 
strong planetary wave driving and weakening of both polar 
vortices during the extended boreal winter season (e.g., 
Tables 2 and 3). 

[ 40 ] The Arctic stratospheric responses to WPEN were 
qualitatively similar under present-day and future climate 
conditions. The weaker magnitude of the late winter Arctic 
response reflected the relatively weaker midwinter forcing 
by SSTs in the Nino 3.4 region in the modeled late 21st 
century SSTs as compared with the present-day SSTs 
(Figures 2 and 12a). Results of an idealized WPEN simu- 
lation suggested that SSTs outside of the tropical Pacific 
have a minor influence on the magnitude of the Arctic strato- 
spheric response. 

[ 41 ] The nonresponse of the Antarctic stratosphere to 
WPEN events, in a future climate, reflected simulated future 
changes in the tropospheric teleconnections. The present-day 
SH response to WPEN was characterized by a strong, pole- 
ward (PSA- 1 -like) wavetrain in the upper troposphere and 
an enhancement of tropospheric eddy heat flux entering the 


polar stratosphere. In the future climate simulations, the pole- 
ward wave response weakened and a zonal wave response 
(PSA-2-like) appeared in austral spring (Figures 4 and 8). 
In this modeling framework, it was not possible to completely 
separate the impacts of future SSTs from those of future 
climate conditions, since the climate conditions feedback to 
the ocean surface and influence the resulting SST fields. 
Nevertheless, sensitivity simulations, in which the prescribed 
SSTs and radiatively active gases were varied in separate 
simulations, provided insight into the qualitative impacts 
of the modeled SST fields and climate forcings. These sim- 
ulations revealed that the simulated future South Pacific 
SSTs, rather than the increases in GHGs and decreases in 
ODSs, caused the change in WPEN teleconnection pattern 
(Figure 12b). Weakening of the poleward wavetrain was con- 
sistent with a marked weakening of the convective response 
to WPEN in the southeastern South Pacific Convergence 
Zone (Figure 10). Similar changes in the upper tropospheric 
teleconnections were found in simulations forced by a differ- 
ent set of modeled 21st century SSTs. In isolation, Antarctic 
ozone recovery had a negligible impact on the polar strato- 
spheric response to WPEN. 

[ 42 ] The appearance of the PSA-2-like pattern and positive 
geopotential height response at 30°S, 120°E, when future 
SSTs were prescribed, reflected strong heating over the 
Australian continent (Figures 5 and 11). While it is beyond 
the scope of this paper to determine the cause of this 
Australian heating or how likely such a heating response to 
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Figure 12. Scatter plots showing the relationship between 
(a) DJF seasonal mean SST anomaly in the Nino 3.4 region 
and the February /March Arctic polar cap temperature 
response at 50hPa and (b) the SON seasonal mean SST 
anomaly in the south central Pacific (180-220°E, 5-15°S) 
and the November/December Antarctic polar cap tempera- 
ture response at 50hPa. Each symbol represents the mean 
difference between WPEN and NTRL simulations: Black cir- 
cles indicate simulations forced by observed SSTs (P, FCLI, 
and A). Grey circles indicate 1DEAL W -P N differences. White 
circles indicate the simulations forced by CCSM3 SST fields 
(FSST, F, and C). White squares indicate B W -B N differences. 
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WPEN is to occur in a future climate, Watterson [2012] 
found that future Australian warm climate regimes are related 
to relatively wanner SSTs in the central and western equatorial 
Pacific and eastern Indian oceans. Furthermore, the PSA-2 
pattern is strongly correlated with SSTs in the Indian Ocean 
and surrounding Australia [Mo and Paegle, 2001]. 

[ 43 ] The above results motivate further study with a 
coupled ocean-atmosphere CCM. Planned simulations with 
a new ocean-atmosphere formulation of the GEOSCCM will 
investigate the Australian climate response to changes in 
nearby SSTs, considering atmosphere-land-ocean feedbacks, 
as well as the stratospheric response to various types of 
ENSO events that are generated interactively by the model. 
SST biases, particularly in the Nino 3.4 and south-central 
Pacific regions, will need to be considered in interpreting 
the modeled atmospheric response to ENSO. 

[44] Acknowledgments. The authors thank NASA's MAP and ACMAP 
programs for funding, and three anonymous reviewers for their helpful feedback. 
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